Quantum Anomalous Hall Effect with Cold Atoms Trapped in a Square Lattice 
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We propose an experimental scheme to realize the quantum anomalous Hall effect in an anisotropic 
square optical lattice which can be generated from available experimental set-ups of double-well 
lattices with minor modifications. A periodic gauge potential induced by atom-light interaction is 
introduced to give a Peierls phase for the nearest-neighbor site hopping. The quantized anomalous 
Hall conductivity is investigated by calculating the Chern number as well as the chiral gapless 
edge states of our system. Furthermore, we show in detail the feasability for its experimental 
detection through light Bragg scattering of the edge and bulk; states with which one can determine 
the topological phase transition from usual insulating phase to quantum anomalous Hall phase. 
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INTRODUCTION 

Twenty years ago, Haldane proposed a toy model in 
the honeycomb lattice to illustrate the quantum anoma- 
lous Hah effect (QAHE) Q, in which a complex second- 
nearest-neighbor hopping term ^26**^ drives the system 
into the topologically insulating state. Different from the 
conventional quantum Hall effect (QHE) Landau lev- 
els are not necessary for QAHE, while in both the QHE 
and QAHE systems, the time-reversal symmetry (TRS) 
is broken. The quantized Hall conductivity can be ex- 
plained with Laughlin's gauge invariance argument [1] 
and by Halperin's edge state picture Q , which has a deep 
topological reason as the first Chern class of U(l) princi- 
pal fiber bundle on a torus [5|] . Considering the topologi- 
cal nontriviality and the absence of magnetic field of the 
QAHE, realizing experimentally this new state of matter 
in its cleanest form is of fundamental importance in the 
study of new materials such as topological insulators. 

Unfortunately, Haldane's model cannot be realized in 
the recently discovered graphene system, since the re- 
quired staggered magnetic flux in the model is extremely 
hard to achieve. A recent proposal predicts the QAHE in 
the Hgi_2;Mna;Te quantum wells \&\ by doping Mn atoms 
in the quantum spin Hall system of the HgTe quantum 
weU to break TRS 0-i]. QAHE is reachable within a 
time range much smaller than the relaxation time of Mn 
spin polarization which is about 10^'^s, while so far the 
experimental study of this effect is not available. On the 
other hand, the technology of ultracold atoms in optical 
lattices allows for a controllable fashion unique access 
to the study of condensed matter physics. An artificial 
version of the staggered magnetic field (Berry curvature) 
with hexagonal symmetry is considered to obtain Hal- 
dane's model for cold atoms trapped in a honeycomb 
optical lattice [13] • While a periodic Berry curvature can 
be readily obtained by couplin g at omic internal states to 
standing waves of laser fields 11 -HI, the experimental 
realization of the staggered magnetic field with hexag- 



onal symmetry remains a challenge. Furthermore, Wu 
shows QAHE can be reached with the p-orbital band in 
the honeycomb optical lattice by applying a technique 
developed in S. Chu's group [17| to rotate each lattice 
site around its own center [l6j. 

In this work, we propose a distinct realization of QAHE 
in a two-dimensional (2D) anisotropic square optical lat- 
tice, which can be realized based on the double- well ex- 
periments performed at NIST [3], superposed with a 
periodic gauge potential which is also experimentally ac- 
cessible. The experimental detection of quantum anoma- 
lous Hall (QAH) states is also proposed and investigated 
in detail through light Bragg scattering. 



THE MODEL 

We consider an anisotropic 2D square optical lattice 
depicted in Fig. 1 (a) filled with fermions (e.g. ^Li, 
'"^K), whose optical potential is expressed as Viatt{i^) = 
-Vb(cos2 kox + cos^ koy) - Vq cos2[^(a; + y) + f ]. This 
potential can be generated from the available experimen- 
tal set-up of^the double- well lattice illustrated in the Fig 
lb of Ref. lLSI by placing an additional polarizing beam 
splitter (PBS) before the mirror reflector M3 and suit- 
ably adjusting the phase difference among different opti- 
cal paths. The first term of the potential is from the light 
component with the in-plane (x-y) polarization which is 
deflected by the PBS and then reflected back by M3, 
while the second term is from the light component with 
the out-of-plane (z)-polarization which passes the PBS 
without reflection. No interference exists between these 
two components. This optical potential has a structure of 
two sublattices A and B. The potential minimum Va at 
site A is higher than Vb at site B as Va — Vb — Vq. 
The anisotropic potential at site A has different fre- 
quencies along the directions of ei.2 — '^{^x i iy) as 

= (Vbfco)^^^/"^^^^ and = V^ujf, and those at site 



B are — ^/Sujf and 



= V^ujf, respectively. 
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The local orbital is described by 2D harmonic oscillator 
eigenstate V'nin, "^ith the eigenvalues 



ni,n2 

Jil ,712 



1, 
1, 



« + -)fi^f + K + ^)f^2 - Vo. (1) 



Below we shall consider Vq taking the value of Vb = 
(SVS - l)'^Er/2 - AM with M satisfying M < Er and 
= h^kQ/2m the recoil energy. In this case, the s- 
orbital at the i?-sites is the lowest one, while the s-orbital 
at the A-sites is nearly degenerate with the p-orbital at 
the B-sites along the ei direction with the energy dif- 
ference of Ef^ - E^^ « 2M. Specifically, if M = we 
have E^Q ~ E^^ ^ Q. For convenience, we denote such 
two nearly degenerate states by V'a — i^oo and ipb — V'foj 
respectively, which consist of a pseudospin-1/2 subspace. 
We shall focus on the hybridized bands between V'a and 
"06 intermediated by the intersite hopping, but neglect 
the hybridization between ijja.t and the lowest s-orbital 
in the B sites since its onsite energy is far separated from 
those of tpa b in the case of large trapping frequencies for 
the lattice m. 
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FIG. 1: (Color online) (a) Fermions trapped in a 2D 
anisotropic square optical lattice. Due to the diflerent on- 
site trapping frequencies, the square lattice is divided into 
two sublattices A and B; (b) The local orbitals are in the ei 
and 62 directions. 

To break TR symmetry, we introduce a periodic 
adiabatic gauge potential in the simple form A(r) = 
?iAosin(fco(y — x))ey with Aq a constant, which can be 
generated by coupling atoms to two opposite-travelling 
standing-wave laser beams with Rabi-frequencies ili = 
nosin{ko{y - x)/2 + Tr/A)e-''^°y,n2 = rJo cos(/co(2/ - 
x)/2+7r/4)e'^o2' [IMl. Note the gauge field can lead to 
an additional contribution |Ap/2m without the square 
lattice symmetry [l9j . which, however, will not distort 
the present square lattice. This is because, first of all, 
the term |Ap/2m is zero at all square lattice sites; sec- 
ondly, the contributed potential of this term is along the 
6*2 direction and leads to the same correction to ^2 ' ■ 
Therefore we still have o;^ = o;^ and the properties of 
■00^6 are unchanged. This is a key difference from the 



situation in honeycomb lattice system, where a gauge 
potential with hexagonal symmetry is required to avoid 
the distortion of original honeycomb lattice [13, m. The 
simple form of gauge potential indicates here a feasible 
scheme in the experimental realization. 

The periodic gauge potential gives rise to a 
Peierls phase for the hopping coefficients obtained by 
exp(i(/)L ) = exp(i J^'~^^ A • dr/h), where the integral is 
along the hopping path from the site to site r^-l-L. Tak- 
ing into account the hopping between both the nearest 
and second-nearest neighbor sites, we obtain the Hamil- 
tonian in the tight-binding form: H = Hi + H2 + Hz, 
with Hz = MErJCfc(i-Oa(r.O - CUr^)Cair^)], and 

ri 3 = 1.2 

+ [E E tabe^''^^Ci{ri)Cb{r, + S,) + h.c], 

i_L—a,b J — 1,3 

+ E t^2e^'^'^^Cl(r,)C,{v,+K,)\, (2) 

3=2A ^ 

where C^(ri) is the annihilation operator on site in 
sublattices A (for fi — a) and B (for jj. — b), the vectors 
Si(-S3) = (a,0),S2(-S4) - (0,a),Ri(-R3) = (a, a) 
and R2(— R4) = (—a, a), with a — n/ko the lattice 
constant. It is easy to check that the Peierls phase 
= 't'Si,, = 0' while (^;^;_^ = 00 (or 0^^^ = -0o) 
with 00 = V2Ao/ko when the site belongs to sub- 
lattice A (or sublattice B). From the symmetry of the 
wave functions ^pa,b we know the hopping coefficients sat- 
isfy tal,ta2,tb2,tab > 0]tbl < 0, tal ^ ta2 and hi ^ tb2 

(Fig. 1(b)). Besides, since the hopping constants decay 
exponentially with distance, tab is typically several times 
bigger in magnitude than taj and tbj {j = 1,2). Nev- 
ertheless, such differences will not affect the topological 
phase transition. Finally, noting ipa,b have the same spa- 
tial distribution in the 6*2 direction, we can expect that 

ta2 ~ tb2 or \ta2 - ^62 | < tal - tbl- 

It is convenient to transform the tight-binding Hamil- 
tonian into momentum space, say, let Ca,b(rj) = 
-^J^k^^^ '^^^a.bO^) and we obtain Hamiltonian in the 
matrix form H = J2]^C^k)n{k)C{k) with C(k) = 
((7a (k), C'b(k))-^ and (neglecting the constant terms) 



n(k) = K{k)a^ + Xy{k)ay + Xz{k)az 



(3) 



Here A^; — 2tab sin0o sin(fcj,a), = 2tab{siii{kxa) + 
cos 00 sin(fcj,a)) and Xz — — M — Ao cos(fc2,a) cos(fcya) — 
2ii sin(A:2;a) sin(A;ya) with ti = {tai — tbi + tb2 — 

ta2)/2, Aq = tal - tbl + ta2 - tb2- As long aS 00 ^ TITT, 

the TR symmetry of the system is broken. This Hamil- 
tonian leads to two energy bands with the spectra given 
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by £± = ±-\/^~Af(k). One can check when ^ mr, 
the band gap is opened if at the two independent Dirac 
points K{Vic) + 0, with (0, 0), (0, tt). Therefore 

when M ^ ±Ao the system is gapped in the bulk. 

When the Fermi energy is inside the band gap, the 
longitudinal conductivity is zero. The anomalous Hall 
conductivity (AHC) can be calculated by 



(4) 



where Ci = (47r) ^ J^^^ dk^dkyU ■ dk^n x 9^^,11 with 
n(k) = (Aa;, Ay, Az)/|A(k)|, is the first Chern number, a 
quantized topological invariant defined on the first Bril- 
louin zone (FEZ). Actually one can construct a map- 
ping degree between FEZ torus and spherical surface 
5^, F : X i-> S^, which gives the Chern num- 
ber Ci = m e Z, with m the times the mapping covers 
the S*^ surface. Ey a straightforward calculation we can 
show in the case — Aq < M < Aq and thus the effective 
masses of the Dirac Hamiltonian around two indepen- 
dent Dirac points are opposite in the sign, the Chern 
number Ci ^ +1 when < (/)o < tt and Ci — +1 when 
— TT < 00 < 0. In all other cases we have Ci = 0. The 
quantized AHC can support topological stable gapless 
edge states on the boundaries of the system. To study 
the edge modes, we consider a hard-wall boundary [21| 
along X axis at a: = and x — L. The momentum is 
no longer a good quantum number, and we shall trans- 
form the terms with k^ in the Hamiltonian back to po- 
sition space. For convenience we envisage first the case 
4>o — '^1'^ s-nd obtain that (H = ^ 'H.(ky^Xi)) 



H ^ 



where 



and 



A4 — 2tab sin kya^x — Muz, 



(5) 



(6) 



A 



i=l,2 



COS kyaaz 



^-^{tai- hi + tb2 - hi ) sin kya<7z . 



(7) 



Note the edge modes are generally exponentially local- 
ized on the boundary [22j. Denoting by ^+(fcy,Xi) and 
'i'^{ky,Xi) the edge states on the boundaries x = and 
X = L, respectively, one can verify that they take the 
following form 



2t2 cos kya — i2ti sin fc^a)) and ^^'"2 = 1 /^j~2' ^^e spinor 
satisfies CxV'''^ = i^'"'': Uk^iy) is the Eloch wave 
function along y axis and Af± the normalization factor. 
The chirality of the edge modes can be found from their 
spectra £^ = ±2tab sinfe^a, respectively associated with 
group velocities vp « zL2atab/fi around Dirac point. Ee- 
sides, the exponential decay of the edge states on bound- 
aries requires that 1^,^21 < 1 ^1.2 > 1 UM- 



Dirac point ky = 0, one can check such inequalities lead 
to "Ao < M < Aq, which is consistent with the con- 
dition for nonzero Chern numbers obtained before. The 
Fig. [2] depicts the energy spectra in different situations. 



E(k)2 




E(k)2 



E(k)2 




(ero^-^l^*, (8) 



FIG. 2: (Color online) Gapped bulk states (black lines) and 
gapless edge states (blue and red lines) on the boundaries 
a; = ((a)(c)) and x — L ((b)(d)). Parameters in (a-b) are 
M = 0,(^0 = Tv/2, and in (c-d) are M = 0, (^0 = -7r/2. The 
chirality of edge states indicates the Chern number Ci = -1-1 
for (po —71/2 and Ci — —1 for (j)o — — 7r/2. 



DETECTION OF THE TOPOLOGICAL PHASE 
TRANSITION 

Next we proceed to study the detection of the edge and 
bulk states with light Eragg scattering, with which one 
can detect the topological phase transition. In the Eragg 
spectroscopy, we shine two lasers on the lattice system, 
with the wave vectors ki, k2 and frequencies wi,2 = A:i,2C, 
respectively (Fig. 3(a)). Note only the momentum ky is 
still a good quantum number, we let q = ki — k2 = qiy 
and denote uj = uii — uj^- The atom- light interacting 
Hamiltonian reads 

H,r^t= Yl ^e''i"'C'liky + q)C,iky) + h.c., (9) 



where the complex variables ^■'"^ obtained by 2 = 
[M ± J A/P + A{tab - iii sin kyu)'^ - Ag cos^ kyo] {2hb - 



where is an effective Rabi-frequency of the two-photon 
process in light Eragg scattering, and the indices /z, v 
may represent the edge or bulk states. The light Bragg 
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scattering directly measures the dynamical structure 
factor (23i]: 



with initial (final) atomic state before (af- 

ter) scattering and f{E) the Fermi distribution function. 
For the topological insulating phase, e.g. when (/)o — 7r/2 
and M = 0, since the edge states are localized on the 
boundaries, we may consider two basic situations for the 
Bragg scattering, say, first we shine the two lasers on one 
boundary (on a: = or a; = L) of the system; secondly 
we shine them on the whole lattice system including both 
boundaries. For the former case, only edge states on 
the boundary shined with lasers can be scattered. When 
Lo < Aq, the initial edge states below Fermi energy will 
be scattered to edge states above Fermi energy, while for 
uj > Aq, part of the initial edge states can be pumped to 
upper band bulk states after scattering. Note that the 
edge state is an exponential decaying function in the x 
direction, with decaying property dependent on the mo- 
mentum ky-^ , while the bulk states are standing waves 
along X axis. Therefore, the scattering process with an 
edge state pumped to bulk states actually includes many 
channels characterized by different values of of the fi- 
nal bulk states, and the effective Rabi-frequency for such 
scattering processes is generally a function of ky^ , kx2 , de- 
noted by ^{ky^, kx2)- Nevertheless, in the practical case, 
we require uj is slightly above Ag, and only the states 
with momenta near zero need to be considered. In this 
way we can expect does not change considerably from 
CIq = n(0, 0), and can be expanded around this value. 
Specifically, in continuum limit one finds the exponential 
decay property of the edge state (localized on x = 0) 
given in Eq. ([8]) can be approximated as ^+ ^ g-gox 
with qo = 2'tt' b ' ^^^h which we obtain the effective Rabi- 
frequency 17 (fca;, fcj^) w fio/(go + ^x)- Bearing these results 
in mind, we can verify for \q\ > qo, the dynamical struc- 
ture takes the following general form 

^(q,^) « {\q\ + kyg)n^SiujTVFq) 
Aq^I^ l + 3M}/{2q()atab) 



,(/) 



ao 
TT 

2+' 



full / {qoatabW 

q-qo 



'fell 



e(c^-c^,),(ii) 

where u = lo — Ao/2?i — 2\q\atab/h,ujc = (Ao/2 + 
2|g|atab)/?i, a — 2Ag^t^ja^, and the step function 
6(x) = 1(0) for X > (x < 0). For q < qo, the S{q,uj) is 
obtained in the same form, only with lUc changed to be 
= (Ao -I- 2a'^tlf,q'^/Ao)/h in Eq. (HU). The first term 
in 5(q, uj) is contributed by the scattering processes with 
edge states pumped to edge states on the boundary a; = 



(for "-") or a; = L (for "+"). The peaks at w = ±qvF 
obtained by this term reflects the chirality of the edge 
modes (Fig. 3(b-c), the blue solid line). From the sec- 
ond term in S'(q, w) we see the scattering processes with 
initial edge states pumped to the upper band bulk states 
approximately give rise to a finite contribution (Fig. 3(b- 
c) , the red dashed line) . For the latter situation the lasers 
are shined on the whole lattice, the cross scattering pro- 
cess with edge states on one boundary scattered to edge 
states on another boundary will happen. We can show 
that the cross scattering leads to an additional contri- 
bution -^Q{uj T vpq) to the dynamical structure (Fig. 
3(c), the blue dotted line). 




x=0 e 



S(q,B) 



(c) 



S(q,co) 





(o=±c|Vp 



FIG. 3: (Color online) Schematic of light Bragg Scattering (a) 
and dynamical structure for the cases that the system is in 
the topological phase (b-c), and in the insulating phase (d). 



Finally, when \M\ > Aq, the Chern number be- 
comes zero and in this case no edge state can sur- 
vive on the boundaries. Specifically, for the case M = 
— Aq (and (j)o = n/2), the bulk gap is closed (A = 
at the Fermi point k = 0) and the particles can be 
described as massless Fermions [24], while for M = 
— 2Ao, the bulk gap is given by A = 2Ao. We re- 
spectively obtain the dynamical structure for the two 



cases S{q,uj) 



S{q,io) 



e{uj 



and 



e(w 



J, with cjci — 2atabq/fi and 



= (2Ao + q'^a'^tlf^/ Ao)/h. The results are plotted 
with blue solid line (for A = 0) and red dashed line (for 
A = 2Ao) in Fig. 3(d). Based on the results of light 
Bragg scattering in different situations, we clearly see 
the Bragg spectroscopy provides a direct way to observe 
the edge states and bulk states. 
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CONCLUSION 

In conclusion, we proposed a novel scheme to real- 
ize the quantum anomalous Hall effect in an anisotropic 
square optical lattice based on the experiment set-up of 
the double- well lattice at NIST. We have studied in detail 
the experimental detection of the edge and bulk states 
through light Bragg scattering, with which one can de- 
termine the topological phase transition from usual insu- 
lating phase to quantum anomalous Hall phase. 
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